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The residual part of the realistic forces — obtained after extracting the monopole terms re- 
sponsible for bulk properties — is strongly dominated by pairing and quadrupole interactions, with 
important or ■ or, octupole and hexadecapole contributions. Their forms differ from the schematic 
ones through normalizations that lead to a universal A~ x ^ scaling for all multipoles. Coupling 
strengths and effective charges are calculated and shown to agree with the observed values. 

21.60. Cs, 21.60.Ev, 21.30.xy 

In the preceeding paper the nuclear Hamiltonian was shown to separate into an "unperturbed" monopole field 
Q\ . "Hm ~ that demands phenomenological treatment - and a "residual" multipole contribution, TLm, that is well given by 
Q\ • existing realistic interactions In principle the task of fully characterizing TL becomes a pure monopole one, but 

t-H \ for this characterization to be of practical use it is necessary to understand the structure of TLm, which is our present 
task. 

It is through empirical evidence that we know that it must be possible to describe nuclei by solving the Schrodinger 
equation: the systems appear to be rather dilute - in the sense that the residual interaction is weak enough to be 
treated in perturbation theory - except for few (sub) shells in the vicinity of the Fermi level. 

The quality of the interactions can be tested convincingly only in cases where exact diagonalizations are feasible 
in spaces large enough (typically one major shell) to ensure that what is left out can be treated perturbatively. The 
number of possible tests - though limited to vector spaces that do not exceed dimensionalities O(10 7 ) - is significant, 
and it is certainly no accident that the realistic TLm passes them all rather well. Therefore, the time may have come 
to trust the realistic forces rather than test them. 

There are many advantages in doing so, but one is special. As new regions open to scrutiny - through Monte 
' Carlo techniques Q and further improvements in shell model technology || - the nature of the problems changes: 
Dimensionalities grow exponentially with the size of the systems, but the behaviour described by the enormous 
*n matrices also becomes simpler through the increasing influence of coherent effects. 

The special advantage in trusting the Hamiltonian is that in adapting existing methods to face this situation, or 
designing new ones, our basic tool is the Hamiltonian itself. 
It is here that a deeper understanding of TLm is called for. 

In section I we are going to show that its structure is as simple as could be expected, because of the strong dominance 
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of pairing and lowest multipole operators. They appear in a normalized form that retains all the simplicity of their 
traditional counterparts but suffers none of their drawbacks. 

In section II, it is explained how to ensure that a given choice of dominant terms is optimally extracted from the 
k^J ' Hamiltonian. 

Section III presents some estimates of the renormalizations due to core polarization effects. 

Section IV is devoted to a comparison of different realistic interactions among themselves, and with the W fit in 
the sd shell §]. 

Appendix A contains results on the matrix defining the diagonal multipole representation. 
Appendix B revisits the pairing plus quadrupole model. 



I. DIAGONAL REPRESENTATIONS OF TL 



Our plan is as follows: a) Introduce TLm in the two standard representations, b) reduce to sums of separable terms 
by diagonalization, c) show thay the few dominant ones are normalized versions of the standard pairing and multipole 
forces, d) show that the normalizations are dictated by the universal A -1 / 3 scaling for all couplings; and finally in e) 
say a few words about monopole the consequences of this result. 

a) The Hamiltonian. We start by borrowing eqs. (1.41 to 1.43) of [Q to write TLm in the normal and mutipole 
representations : 

TL M = W? stu Z+ r -Z tur , (1.1a) 

r<s,t<u,T 
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hm = £ M 1/2 (1 + 6rsf, J + ^" )1/2 ^..row, (i-ib) 



^j(-r +t - 7 - r {; ; ^ Kun. ti..],. 



r 

< stu = E(-) S+t_7_r { II 7 }^taJ7], (1-ld) 

Since H m is defined as containing all the 7 = 00 and 01 terms, Hm is defined by 

^rstu = f° r 7 = 00 and 01 

which explains the absence of one body contractions. 

The normal or W^-form is unique and the ordering of the indeces simply eliminates double counting: the contributions 
in rtsu,rstu,rsut and rust are identical, and it is just as well to keep only the first. 

The multipole or w-form is highly non unique because the terms are not linearly independent and permuting indeces 
leads to different objects. We have chosen the variant in which summations are unrestricted for a reason that will 
become immediately apparent. 

For the calculations we adopt the KLS force , but in section IV the claim that the choice of realistic interaction 
does not matter much. 

b) Separable form. We call Hm the restriction of Hm to a finite set of orbits. Replacing pairs by single indeces 
rs = x, tu = y in (1) and rt = a, su = b in (2), and bringing the matrices W x and f2 b = w ab\/(l + S rs )(l + <5t u )/4, 



to diagonal form through unitary transformations U xk , 



ak ' 



U- l WU = E^W T xy =Y J U T xk Ul k E T k (1.2a) 



« - 7« = e— ►& = £«Zk«Z* e Z' (1.2b) 

k 

and then, 

H m = E E l E U ^r ■ E U yk Zyr, (l-3a) 

fe,r x y 



% = E e lE "a A 7 E < s 2 M 1/2 , (i.3b) 

fc,7 V a b ) 

which we call the E and e representations. Note here the explanation of the unrestricted ordering of the orbital 
indeces: it guarantees that in the j a b matrices, a and b belong to the same set. In Appendix A it is explained what 
happens when they do not (asymmetric factorization). 

c) Dominant terms. We have calculated the eigensolutions in (4) using KLS for spaces of one and two major 
oscillator shells. The density of eigenvalues (their number in a given interval) in the E representation is shown in fig. 1 
for a typical two-shell case. It is skewed, with a tail at negative energies which is what we expect from an attractive 
interaction. 

The e eigenvalues have a number of simple properties demonstrated in Appendix A: their mean value always 
vanishes, their width is y/T/8 of that of the E distibution, and they are twice as numerous. In fig. 2 we find that 
they are very symmetrically distributed around a narrow central group, but few of them are neatly detached. The 
strongest have 7^ = 1 _ 0, 1 + 1, 2 + 0, 3 _ 0,4 + 0. If the corresponding states are eliminated from H in (4b) and the 
associated H in (4a) is recalculated, the E distribution becomes quite symmetric. Details will be given in section II, 
and here we only note that the residual skewness is entirely accounted for by the V — l + 0, + l and 2 + peaks, of 
which the first remains strong at — 7MeV. 
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This result is most telling because from the work of Mon and French Q we know that a symmetric E distribution 
will lead to spectra in the n-particle systems that are identical to those of a random matrix. Therefore, we have found 
that - with the exception of three T peaks - the very few dominant states in the e-distribution are responsible for 
deviations from random behaviour in Hm- Positively stated, these states are at the origin of collective properties. 

If the diagonalizations are restricted to one major shell, negative parity peaks are absent, but for the positive parity 
ones the results are practically identical to those of figs. 1 and 2, except that the energies are halved, a striking feature 
whose significance will become clear soon. 

In the list of important contributions whose structure we analyze we include the T = 10 and 01 terms, and the six 
strongest 7 ones 

Their eigenstates (i.e. the factors in eqs.(4) with k = 1) will be compared with standard pairing and multipole 
operators. To fix ideas we write the form these eigenstates should take in the case of perfect pairing (T = 01) and 
quadrupole forces (7 = 20) acting in one shell of principal quantum number p. To compare with the result of a 
diagonalization, the operators must be normalized: 

Pp = Poi P = E Z -oi^ /2 /^p /2 , % = 0.655 A 2 / 3 (1.4) 
q p = Mf = J2 SlUrs/Afp, Np = 0.085 Af/ p 3 (1.5) 

rsGp 

where 

• n r = j r + 1/2, q rs = < r||r 2 r 2 ||s > j sjl 

• A m p is the the total number of particles at midshell p (p^ = p(p — 1), remember) 

A mp = 2^(p' + 2)( 2 ' + (p + 2)( 2 ' = 
p'<p 

i((p + l)(p + 2)(2p + 3) « H(p + 3 /2) 3 . 

• the norms f2 p and Af p , are then 

n p = J2^r = (p + 2) (2) /2 « \ (3A mp /2) 2/3 , 

r 

Ml - Sg 2 s - 5(p + 3/2) 4 /32tt = J_( 3j4mp /2) 4 / 3 . 

OZ7T 

For the calculation of Af£ we have used the matrix elements listed in (All-15) of [1]. 

For the other strong multipoles the choice of operators is evident and for T = 10 the simple idea is that pairing in 
LS coupling should produce a good candidate. Labelling the orbits by their £ quantum numbers, we have two pairing 
terms 

Ptr = (EM)~ 1/2 EM 1/2 ^W ST = 01, 10, 

1 i 

which in jj coupling become = f2~ \/Q^Z^ r01 

r 

P io = ^ 1/2 Y,W 2 {{U)Q{1/21/2)1 \(jj')l}Z+ r , w , (1.6) 

where we have recovered the usual r = j£ label, and used a self evident notation for the LS to jj transformation. In 
table I these operators are compared with the results of the diagonalization. It is apparent that P p oi accounts very 
well for (U~rr)p- For P p iq vs (U^) p the agrement is not so excellent, but still good. The overlaps are found under 
< U p \P >) in table II, which also contains the corresponding values of < u p \M > for the lowest multipole operators 
M. The agreement is again excellent except for the a case, for which it is only fair. Note that the form of these 
operators is given in the Appendix of [1]. 
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TABLE I. Eigenvectors and energies calculated in the pf(p = 3), sdg(p — 4), and pf + sdg(3 + 4) spaces compared with the 
normalized pairing operators Ao and Poi 





Poi 


(Crs)3,4 


^2(^)3+4 


Ao 






77 


-.63 


-.65 


-.66 


.41 


.31 


.31 


75 








-.68 


-.74 


-.74 


33 


-.45 


-.38 


-.41 


.33 


.26 


.27 


35 








0. 


-.30 


-.25 


31 








-.42 


-.43 


-.44 


55 


-.55 


-.58 


-.59 


-.27 


-.12 


-.13 


11 


-.32 


-.29 


-.31 


-.04 


.07 


.04 


-Et(3) 




-2.95 






-4.59 




99 


.58 


.61 


.64 


-.37 


-.26 


-.25 


97 








.63 


.67 


.72 


77 


.52 


.55 


.57 


.26 


.12 


.11 


75 








0. 


-.27 


-.28 


55 


.45 


.41 


.34 


-.31 


-.22 


-.20 


53 








.46 


.50 


.48 


33 


.37 


.35 


.28 


.16 


.03 


.01 


31 











-.21 


-.18 


11 


.26 


.17 


.13 


-.26 


-.22 


-.20 


-Et(4) 




-2.65 






-4.78 




E r(3+4) 

2 






-2.76 






-5.06 



TABLE II. Energies and eigenstates of the dominant terms (7 = 21 added for illustrative proposes). See text. 



7 


e 3 


p 1 

e 4 


p 1 

e 3+4 


M 


< u 3 \M > 


< u 4 \M > 


< M3K1 > 


< 114 \u' 2 > 


a 2 


11 


1.77 


2.01 


3.90 


GT 


.992 


.994 


.999 


1.000 


.94 


20 


-1.97 


-2.14 


-3.88 


r 2 Y 2 


.996 


.997 


1.000 


1.000 


.95 


10 


-1.02 


-0.97 


-1.96 


a 


.880 


.863 


.997 


.994 


1.04 


21 


-0.75 


-0.85 


-1.60 


r 2 Y 2 T 


.991 


.998 


.999 


.997 


.94 


r 


El 


El 


E3+4, 


P 


< U 3 \P > 


< Ua\P > 


< Uz\U[ > 


< Ui\U^ > 


a 2 


01 


-2.95 


-2.65 


-5.51 


Poi 


.992 


.998 


1.000 


.994 


1.048 


10 


-4.59 


-4.78 


-10.12 


P10 


.928 


.910 


.998 


.997 


.991 
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Interesting as these results might be, the truly remarkable ones come when we diagonalize in two major shells. Let 
us go back to table I, and note that the eigenstates can always be written as 



C/3+4 = (aC/g + PUfi/yfi a 2 +(3 2 = 2, (1.7) 

j' v can - in principle - be any unit vector, but in fact it is almost identical to U p . This is always the case as 
table II shows: < Vi\U[ > and < u^u^ > are strikingly close to 1 with no exceptions while a 2 is quite close to 1. 
Therefore, for any normalized pairing or multipole operator we have that 

If O p and Op+i are eigenstates for shell p and p+1 separately, then (O p + Op+i) is very much the eigenstate for 
the space of the two shells. The eigenvalues are very close in the three cases. 

Note that we have chosen a normalization of 2 for the two shell eigenstates so as to halve its eigenstate. 

Before we examine the consequences of this result, we mention a few facts about the other contributions. 

In table II we have added the 7 = 21 case as a reminder that isovector multipoles are always present. Their strength 
is between 30 and 40% of that of the isoscalar terms (precise numbers are given in table III in section II) and they 
have identical structure. 

The dominant negative parity contributions are 

7 = 10 at 4.59 MeV, < rY^u 10 >= 0.994 

7 = 30 at 2.69 MeV, < r 3 Y 3 \u 3a >= 0.986 

The first is a center of mass operator. Its presence simply reflects the translation invariance of the interaction. Its 
7 = 11 counterpart, associated to the giant dipole resonance (GDR), comes at 1.81 MeV. The other strong term is 
responsible for octupole collectivity. (In deciding whether a given multipole is attractive or repulsive it should be 
remember that (M 7 M 7 ) = (-) 7 [7]~ 1/2 (M 7 • M 7 )). 

d) Universal scaling. We have now the necessary elements to construct a schematic but accurate collective 
Hamiltonian. From 101 we know that 



Wl y {w) S —W£ y (uo) (1.8) 



CO 

and therefore the eigenvalues in table II must scale in the same way. Setting a 2 = 1 in eq.(1.7) for simplicity, the 
normalized pairing and quadrupole forces become 

Hp = -^- o \E°i\(P+ + P+ +1 ) ■ (P P + P P+1 ) (1.9) 

Hq = -^\e 2Q \(q p + q p+1 ) ■ (q p + q p+ i). (1.10) 

IIUJQ 

which we take as representative of the "collective" Hamiltonian because of their known coherence. For the other 
strong terms the expressions are strictly similar, and all arguments concerning pairing and quadrupole expressions 
apply to them. Since ej = ej — eJ +2 /^ (same for E 1 ), the coupling constants could be taken to be independent of 
the space chosen; which may be any of the shells or the two together. 

The term normalized applies to the one shell operators. For two or more shells it is more convenient not to normalize 
their sum, since it is in this form that the couplings are constant. 



Equations (1.9) and (1.1C ) call for a generalization to an arbitrary number of shells and it does not take much 
imagination to discover what it should look like. However it is interesting to understand the origin of this very welcome 
simplicity. Let us start from the obvious it is always possible to write the lowest state associated to multipole 7 as 

To say something rigorous about the radial form of O 7 does not seem trivial, but granted that it is close to the 
multipole operator in shell p, then A/" 7 is fixed by a scaling argument. 

In leading order the expectation value of a Hamiltonian must go as the number of particles in the system. Therefore, 
in 7i, the leading monopole terms must go as 0(A), and TLm acting on one shell must go at most as 0(D = 0(A 2 / 3 ), 
where D is the degeneracy of the Fermi shell. Since there are p possible contributors, each individual multipole term 
must go as 0(A 1 ^ 3 ). 

It is precisely what the normalized operators ensure, given the universal scaling provided by the huj factor. 

In our particular cases: n particles in shell p can produce an energy proportional to nD for the pairing, and n 2 D for 
the quadrupole forces, which the normalizations reduce to n and n 2 / D. Since n = O(D) = 0(A 2 ' 3 ), when multiplied 
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by Two = 0(A~ 1 / 3 ) both contributions become 0(A 1 ^ 3 ). If the n particles are promoted to some higher shell with 
p = p + M, there is no pairing gain, a slight quadrupole loss and a monopole loss O(Mlu) that remains 0(A 1 ^ 3 ) as 
long as M is not too large. 

For the conventional pairing and quadrupole forces the energies in one shell would be nD/A and n 2 D/A 5 / 3 respec- 
tively, as their scalings are chosen to yield the correct 0(A 1 / 3 ) magnitude. But now when the particles are promoted 
D — ► (p + M) 2 , and not only there is a gain: it becomes larger than the monopole loss that is only linear in M and the 
system collapses. If the forces are restricted to act in finite spaces, to obtain sensible results the coupling constants 
must be reduced as the space is increased. (The problem will be discussed in Appendix B). 

The normalization is defined only to leading order 0(p k ), and we can say nothing about the 0(p k ^ 1 ) terms, 
responsible for a 2 ^ 1, and the slight differences we have neglected in the e and E couplings, but this is the good 
reason to neglect them. 

What the argument does not explain is why the pairing and multipole operators in one shell resemble their ideal 
counterparts as closely as we have found so far. There is a hint, however, in that they must be the ones capable of 
producing coherence. More on this on the next subsection, and in section III, it will be shown that the q' p that couples 
shell p to 2Uuj jumps, remains indeed close to q p . 

We have found therefore, that we can recover the geometrical simplicity of the pairing plus quadrupole model 
without its fundamental flaw: the space dependence of the coupling constants. The model has an enormous historical 
interest, and it is very instructive to show how far we can go in justifying it (see Appendix B). 

e) The monopole hint. If Hm, is as good as we have argued why not trust the information it can provide about 
H m 7 It is quite possible that, rather than wrong, it is only insufficient and its study is most interesting. 

As befits the leading term in a multipole expansion, the monopole one is the strongest: in fig. 2 it would come 
at —10 MeV. Unsurprisingly it has the form X) n pl yDpi which is what we expect of normalized operators, but it 
should be mentioned that there are several monopole candidates one can think of: 1, which is the h operator in second 
quantization; but also r 2 for instance. When normalized they give the same result. 

The remarkable thing about this form is that it provides the answer to the problem raised at the end of the 
preceeding paper. Which is that fin, suffers from the defect of the conventional separable forces in that it must be 
associated with a coupling constant that is space dependent. There is no collapse now because n is a conserved 
quantity, and an A^ 1 scaling ensures the asymptotically correct behaviour and the coupling tends to a constant. 
Contrary to the pairing and quadrupole cases there can be no energetic gain in promoting particles to higher orbits 
but there is no loss either, and this is a subtle form of collapse because a good Hamiltonian must ensure the existence 
of a Fermi level, i.e., it must force the particles to occupy the lowest orbits. The normalized monopole operator does it 
by producing a discontinuity at each shell closure. Therefore, it is not only responsible for the bulk energy of nuclear 
matter, but it also takes care of the major shell effects. 

This operator suggests the starting point in the construction of H m . For a preliminary attempt see p5| , where a 
mass formula - of rather high precision by present standards - is derived. 



II. CHOICE AND COMPLETE EXTRACTION OF H c 



The choice of He- The results so far invite a separation: 

H = H m + He + H r 

where He is the collective or coherent part, while Hp is the rest. To define them with some precision , we shall rely 
on the result of Mon and French || that Hp could be viewed as random , as soon as its E distribution becomes 
symmetric. 

The distributions will be characterized by their moments: 

2 

m k = D ( D _^ J2 E r> l3 = m 3 /m 3 2 /2 

and the vanishing of the skewness 73 will be seen to be sufficient to ensure symmetry. We use m-i = a 2 . 

Let us then define some cut off e, eliminate from H in the e representation those states with |e 7 | > e, and decrease 
e until the E distribution becomes symmetric. The result of the operation is shown in figs. 3,4 for e —2 and 1.3 MeV 
respectively. Labelling the original distribution in fig. 1 as (e = 00), we find the following moments (a 2 in MeV) 



G 



TABLE III. The sdg + pf states in the /-representation (|e 7 | > 1.3) and their pf and sdg counterparts. S n , A* = symmetry 
type and parity, fr signals the states singled out in fig.2, and their one shell counterparts. [ is for states with e 7 | < 1.3MeV 
that are likely to have a clear multipole character. Parenthesis indicate that the assignment is unchecked but given as plausible. 



At 


sym 


type 


Spfsdg 


e Pf 




00 


A~ 




1.88 










10 


S~ 


riYi 


4.59 


fr 










S+ 


a 


-1.96 


tr 


-1.02 


-0.97 


fr 




s~ 




-1.53 












s+ 


I 


1.44 




0.66 


0.80 






s~ 




1.41 










11 


5+ 
A~ 
S~ 


OT 

riFiT 


3.90 
-1.83 
1.81 


fr 


1.77 


-2.01 


fr 


20 


s+ 

S+ 


r 2 Y 2 


-3.88 
1.31 


fr 


-1.97 
0.64 


-2.14 
0.75 


fr 


21 


s+ 
s~ 

A 


r 2 Y 2 r 
M2?? 


-1.60 
-1.55 
1.46 




-0.75 

0.64 


-0.85 
0.76 




30 


s~ 


r 3 Y 3 


2.69 


fr 








31 


s~ 


r 3 Y 3 T 


1.14 


1 








40 


s+ 


(r 4 Y 4 ) 
(r 4 Y 4 r) 


-2.11 


fr 


-1.12 


-1.24 


fr 


41 


S+ 


-0.91 


l 








50 


s~ 


(r 5 Y 5 ) 


1.75 










51 


s~ 


(r 5 Y 5 r) 


0.78 


l 








60 


s+ 


(r 6 Y 6 ) 


-1.26 


l 


-0.73 


-0.82 
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• e = oo er 2 = .99, 73 = -2.22. 

Since the lowest state T = l + is at —10.12 MeV its contribution to 73 is by far the largest but still only —0.61. 
It means that many states in the tail must contribute to 7713. 

• e = 2.0 a 2 = 0.60, 73 = -0.79. 

Five peaks have been excluded 7 = 1^0, 1+1, 2+0, 3~0 and 4+0 and now the T = 1+0 state at -7.79 MeV 
accounts for most of 73 with a contribution of —0.60, which when added to that of the next two states T = 2+0 
(—4.29 MeV) and 0+1 (—3.82 MeV), at 73 = —0.86, exceeds the full value. As anticipated, the very large e states 
are responsible for most of the tail of the original distribution. Although distorsions are still apparent, they may 
be interpreted as fluctuations and we have a good model for He = the five 7 states and the T = 1+0, + l ones. 

Since it is not only the size, but also the ability to generate coherence that must characterize the main terms 
we have left out T = 2 + which - as we shall see in section IV - must be counted as a nuisance, rather than a 
bona fide candidate to Ho- 
lt should be noted that the structure of the large T states is little changed in going from e = 00 to e = 2. 

• e = 1.3 a 2 = 0.41, 73 = -0.002 

This cut-off was chosen as a sensible definition of the bulk of the e distribution. The histogramm in fig 4 now 
becomes structureless. The 7 states with |e 7 | > 1.3 are listed in table 3. Some useful information is given, as 
their symmetry type (see eqs. (1.30-32) in [1]), and their multipole nature whenever possible. 

The cut off is reasonable in that the majority of excluded states belong naturally to what is expected from a 
multipole decomposition. However, several peaks have an ill defined status, while others we would include in He have 
missed the cut off (e.g. 3~1 at e = 1.14). Clearly we are in the boundary where the small terms in He overlap the 
large ones in H R . 

Two attitudes compete in choosing He- Include either as much as necessary or as little as possible. Which one is 
given precedence depends on the computational strategy adopted. 

As much as necessary applies to the recently developped shell model Monte Carlo approach The Hamiltonian 
must be written in some e representation, and the authors have chosen an ingenious one in which they introduce the 
Pauli violating terms necessary to cancel all 7 = Al terms. To understand how it can be done refer to eq. (1.47b) 
in the precceding paper [1], and imagine that in V/ S J U the rs and tu states are allowed to be symmetric. Since the 
purpose of the exercise is to reduce the number of terms, it is probably a good idea to proceed as we have done and 
choose a proper He- The problem is that pairing has to be included in the e representation which may be costly, but 
an example given at the end of the section suggests that the cost is probably not too high. 

As little as possible, can be reconciled with as much as necessary by treating first the truly important terms and 
then doing some form of perturbation theory to account for the rest. Here pairing plus quadrupole or even only one 
of them may be sufficient. 

One example may be mentioned. In a forthcoming study fll|| , it is shown that states of 4 neutrons and 4 protons 
moving in the same or contiguous major shells, generate rotational spectra exhibiting systematic backbending. With 
a proper q ■ q the wavefunctions have overlaps of better than 0.95 with those of the realistic interaction H , but the 
backbending has gone. However when the expectation value of the full H wis calculated using the q ■ q eigenstates, the 
spectrum reproduces perfectly the exact one: lowest order perturbation theory may work very well with a reasonably 
good starting vector. 

This example is particularly relevant because the "proper" quadrupole that reproduces best the exact results is 
some 30% larger than e eigenvalue in table 3. It is probable that part of the effect is due to the isovector quadrupole 
that was not included in the calculation. But there is another cause as we see next. 

Complete Extraction 

He and H R in the normal representation, can be thought as vectors, each matrix element W^ stuC and W^ stuR , 
being associated [r] to unit vectors. Then the norm is simply 1112 = cr 2 (assume in all that follows that the matrix 
elements have been all divided by D (D + l)/2) 

a 2 =< (H C + H R ) 2 >-]T« st J 2 [r] 

= E [( W rstuR) 2 + « st uc) 2 + ™l tuR Wl tuC \ [T] 

= H 2 +H 2 +2H C -H R (2.1) 
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It is elementary to show that if Hi and H 2 are the Hamiltonians associated to eigenvectors E\ and E 2 , then 
< H X H 2 >= Hi-H 2 = 0. 

The crux of the matter is that He and Hr are by construction orthogonal in the multipole representation, but 
they need not be in the normal one. We may even invent vectors orthogonal in the former that become parallel 
in the latter. This, because of the non linear independence problem encountered several times in [1]: a multipole 
decomposition can teach us what to separate, but not how to separate. 

What has been done in the preceeding section illustrates quite well the first part of the proposition, and incidentally 
the need to use the / matrix elements in eqs.(1.2), to ensure symmetric factorizations, since any other choice would 
have led to opaque results. Now that we know what to separate, we should do it well. 

The simplest way is to write 

H r + Ho=(h r -%^Hc)+Ho(i ' HC ' HR 



He • He J V He ■ H 

which makes the two terms orthogonal. By squaring we find 



H 



(H C ■ H R f 
{H c ■ H c ) 



Hi + 2H C ■ H R 



{He ■ Hr? 
{H c ■ H c ) 



an since Tg^rg^j can be safely neglected, the orthogonalization amounts to leave Hr almost untouched, and to boost 



He by the cross term. Now we call upon eq. (A6), which says that the total norms in the two representations are 
proportional and examine the relative contributions of He and Hr in both cases. 

For e = 2 , from table 3 we can find that the five largest peaks contribute .29 , i.e. some 30% to the total norm 
which we know to be .99, and therefore Hr accounts for the remaining 70%. But we also know that - when calculated 
in the normal representation - H R is not .70, but .60, and therefore the orthogonalization boosts the share of He to 
40%, which implies a \/4/3 = 15% increase in the coupling constants. 

This should be sufficient to have an idea of the effect, but for a quantitatively reliable extraction, each term must 
be treated separately. Let us see how. 

Assume we have selected the candidates to He, transformed to the normal represetation the multipole ones, and 
made them monopolc free. Call them H K . We want to find the linear combination He = ^2C K H K , that maximizes 
the overlap between He and H . It amounts to solve the standard problem (a = rstu) 

E« " E CM.m = (H - E C K H,f 

a K K 

= H ■ H - 2J2C K H K ■ H + 2J2C K C K ,H R ■ H K , - min (2.2) 

and therefore C K is determined by the linear system 

H K -H = Y J H K -H K ,C K , (2.3) 

k' 

The H K vectors can be made orthogonal by diagonalizing the norm matrix H K ■ H K i through transformation T^ u 
Hfj = E T^Hk , H K = E Tk.^H^ ; Hfj, ■ H^ = 8^1 H^ . (2.4) 
Therefore, defining He in terms of the H^, we have 

Hc = J2 C 'h h k ^ C tl = H^-H/Hl. (2.5) 

As we see the analysis in the first section has fufilled its part in pointing to the correct forms to be extracted, but 
the associated coupling constants may need some corrections. If our estimates are correct, they should not exceed 
20%. 

It is worth emphasizing that once the fully extracted He is transformed back to the multipole representation it 
may look quite different from the original, which is not monopole free, and contains Pauli violating terms. However, 
it is very much the original. 

A potentially useful variant of this possibility was found in the sd shell: When the strongest contribution for each 
7 is included and extraction is carried maximizing the overlap, the addition of the pairing term makes no difference. 
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III. CORE POLARIZATION 



To account for the effect of states outside of one or two major shells in which we are prepared to work, we use the 
quasiconfiguration method jl2|,[l3]]. The idea is to separate the full space into a model space containing states \i > 
and an external one made of states \ j >. The \i > and \j > states are now "dressed" through a transformation that 
respects strict orthogonality 



h >= \i > + 



J2 Ai ^> \j>=\j > -J2 Ai ^> <*ij>=°- 



The amplitudes Aij are then defined by 



< i\H\j >= 0. 



Perturbation theory is the natural way to determine Aij, since whatever must be treated exactly is in principle 
contained in the model space, which we take here to be a full shell p. What is nice about the quasiconfiguration 
method, is that it makes it very clear how the dressing of the states is transformed into the dressing of the operators. 

We refer for details to |]l^ , |l3| , from which we borrow the expression for the second order correction, and apply in 
some detail to the pairing, and quadrupole Hamiltonians. 

Pairing. We deal first with the effect of shell p+1 only. £7oi is a one shell energy from table 2, scaled by hoj/hojo, 
and hui = 9MeV. 



< i\Hp\i >=< il^oiM — 



P+ ■ Pn 



it in 



f > 



(3.1a) 



BgiH ^ < i\P+ • P P +i\J >< J\P+ +1 • P P \i' > 



2Huj 



(3.1b) 



=< i\E 01 (u)) 



p+ . p 

r p r p 



Ej 1 (Lu)lP p + -PpPp + i-P+ +ll . / 



2Huj 3 Qr, 



n 



p+i 



I* >= 



(3.1c) 



< i 



1 



p+ . p 

r P r p 

i2 n 



\i > 



(3.1d) 



Step by step we find 

a) the unperturbed energy, to which we add 

b) the second order perturbation. The intermediate states are assumed to be all at energy 2Tilu. Then 

c) invoking closure and recoupling the operators (hence the factor 1/3), 

d) the final result follows by contracting out the the operators in shell p+1 (which gives back a factor 3). 
Repeating the operation to account for states in shell p — 1 leds to exactly the same correction and we conclude 

that the renormalized operator is the original one affected by a modified Eqi 



Eqi 



E 01 1 



\E,i\(cj) \ _ 
r | = Eqi 

noj 



£7oi(H-0.32), 



(3.2) 



Now consider what would happen if we were to generalize to an indefinite number of shells, and compare with the 
result for ordinary pairing: 



E q1 (lo) — ► £7 iM 
\G\ 



1 + 



l^oi I M 

2huj 



1 



1 1 



G 



G 



1 



2hu 



p+1 



+ 



p+2 



+ 



p+3 



(3.3) 
(3.4) 



The normalization of the operators P p transforms a bad divergence into a logarithmic one. Since high lying pairing 
excitations violate badly translation invariance the corresponding matrix elements should be quenched, and this 



10 



is probably sufficient to eliminate the divergence. Plausible as this argument may be, the fact remains that the 
renormalization of the pairing force is a delicate problem that has yet to be solved. Some more on this at the end of 
the section. 

Quadrupole. The problem of multipole renormalizations is very different from the one we have analyzed for 
pairing because the physical processes are different. Where we had two particles jumping to higher shells, or being 
promoted from lower ones, now it is a particle-hole excitation that is produced by the valence particles, and it is in 
the multipole representation that the problem is naturally treated, and the space involved is always finite in second 
order for the operators of interest. 

In particular, the quadrupole renormalization are due to 2h~uj jumps mediated by a term of the form 
' I'lhu) l'iKu ' %)i which must be extracted from Jiu- The primes indicate that the operators may not be 
genuine quadrupoles. The extraction method is a variant of the diagonalizations we have used so far. It is explained 
in Appendix A under the item asymmetric factorization. The space involves all shells from p — to 4, and the matrix 
elements interested in are eliminated (those of the form q' 2hul ■ g^u) Note that there is no cheating here: we could 
have done the same thing in the two shell case by keeping only the cross shell elements if we had wanted only the 
q p +i component. To compare with empirical reults in the sd, p = 2 shell we work at Tilu = 11. The standard one 



shell calculation yields e 20 = — 2. 40 Me V (note that eq.(B4) gives — 2.37MeV). For the cross terms the large matrix 
produced 

fc = -2.80MeV, < q' 2 \q 2 >= 0.97, (3.5) 
< «»iJ<ftn + 9i3 + 924 > /N 2 hu = 0.83 (3.6) 

where q pp ' are the 2tko quadrupole operators that can couple to q 2 . Their sum is normalized to M 2 % u . By normalizing 
each q vv i ope rator separately the overlap at 0.87 is slightly better, but it is the number in ( |3.6| ) we shall need. 



From (3.5) q' 2 ~ q 2 , and this is an important message: it does not matter much to what q 2 is coupled - even Afuj 
jumps if Ti. demands them (not very likely) - but it must be q 2 to ensure that "quadrupole renormalizes quadrupole" 
and that effective charg es a re state independent as shown below. 



Proceeding as in eqs.( |3.l| ) , we have (he stands for Hermitean conjugate of the first term) 

< t[H q \t >= e 20 (uj) < i\q p ■ q p \i' > - (3.7a) 

< ' % + hc )\ J >< J \(Vm u ■ % + hc)\%' >= (3.7b) 

i , \- - 2k 2 (uj) 1 . . _ _ , 
=< i\e 20 (u>)q p ■ q p —— -q 2hu) ■ q 2HaJ q p ■ q p \i >= (3.7c) 

. . k 2 (uj)\ , , , 

e 2 oM - J < i\q P ■ q P \i > (3.7d) 

In the last step only half of the operators act: those creating first an exitation and then destroying it. Hence the 
factor 5/2 instead of 5. We have indulged in the fallacious approximation of treating the three terms in q' 2hul = 
Q P ~2 P + 9p-ip+i + l' PP +2 as commuting with q p , which is true only for the middle one. A correct calculation would 
yield a rank 2 + 3 force for the offending terms. Still the result is correct for the two body contribution as can be 
checked in ]l^| in which a very similar case is fully worked out. The neglect of three body contributions is common, 
but bad, practice. 

The modification of the transition operator q p (not q p ) is calculated along similar lines, 

< i\q\i' >=< i\q p \i' > -^r^ < i\q2hu,(<l2hu ■ % + hc )\i' > ( 3 - 8a ) 

. -1 I •/ . 2fc(u) i Nzhu „ .| - / ,./. zoom 

=< i\q p \i > ^7~5^v^ < l \lpQ2hu ■ q 2 nJ l > ( 3 - 8b ) 

-°- 83 ^^f) < ^ |i ' > (3 - 8c) 
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what has been done is: 

a) use the second order expression. 

b) upon recoupling collect 4 equal terms and the factor 1/5. Then, to contract the normalized 2tuv operators (whose 
ovelap is 0.83), interchange normalizations for q p and q2h^, and write 

c) the final result. 

Inserting the necessary numbers (A/mw/A/J, = 1-97 for p = 2), the transition operator is boosted by a factor (1+0.42), 
which is too small to agree with the empirical value close to 2. The problem is that second order perturbation may 
work very well for the energy but poorly for the transition rates [fl2f . We shall not go into the reasons, and simply 
borrow a nice argument due to Mottelson ( JlS"]). The idea is that the second order perturbation should not be 
understood as affecting only the valence particles but the whole system. Then, calling Q 20 the total quadrupole 
operator ,the estimate Q 2Q = (1 + 0.42)g 20 becomes 



q,u = q zu + 0A2Q zv q20 = (1 _ o.42)-Y°, 

a result equivalent to an RPA resummation. 

The comparison with empirical data is rewarding: ( rf 1 are the effective charges) 

e 2o — , e 20 (i + 0.3) = -3.12 us -3.18 [16] (3.9) 
r/ 2Q = (1 - 0.42)- 1 = 1.76 vs 1.78(3) § (3.10) 
rj 21 = (1 + 0.19)- 1 = 0.8 vs 0.8(1) § (3.11) 

(The isovector effective charge can be calculated knowing k=1.28) 

There is perhaps an element of luck here, but the presence of the overlap and the norms in the last equation fl3.8| ) 
indicates that the diagonal representation is probably incorporating some important effects neglected in previous 
calculations of the effective charges |lj . 



The comparison of the effective pairing numbers in (3.2), with the empirical values would also be satisfactory, but 



fortuitous, because of the logarithmic divergence of second order theory. 

This problem does not seem to have attracted much attention, but there have been numerous controversies in the 
literature on the influence of core polarization in the strict sense - i.e. what we call multipole processes - on pairing 
renormalization (see fl4j| for a recent review) . 

According to our simple ideas "pairing renormalizes pairing" and "multipole rcnormalizcs multipole" and the 
controversies would seem without object, but there is a catch. 

A pairing force can always be written in the multipole representation, the highest multipoles entering with greater 
weights. In the sd shell - the region that generated the controversies -, quadrupole and hexadecapole are already 
high multipolarities, well represented in the pairing decomposition, and their strong polarizabilities will have a non 
negligible effect. This is simply a manifestation of the lack of linear independence we have often encountered. 

In heavier nuclei the important contributions to He will become increasingly linear independent, and they will 
increasingly renormalize themselves, and the question will remain: what is the precise nature of pairing renormaliza- 
tion? 

Whatever the answer to the question, it is clear that the use of diagonal representations simplifies considerably the 
calculation of effective operators by making it possible to concentrate attention on the few important terms that grow 
bigger, and neglect the small ones that have no chance to grow. 



IV. COMPARING INTERACTIONS 



Now we compare several realistic interactions in the sd shell, the only region in which a direct fit to the data, W 
H , leads to more accurate spectra. 

Table 4 shows some off diagonal matrix elements at the beginning of the shell, calculated with the KLS (as used in 
||), KB [ p!8[ , and Bonn B |l!| realistic interactions, or taken from the VF-fit. It is apparent that the realistic values 
are close to one another (especially KB and Bonn B) and not far from W, except in two cases. 

The comparison of individual matrix elements can be misleading, because apparently large discrepancies may be 
of little consequence, and apparently small ones disastrous. Therefore it is better to concentrate on the contributions 
that can make a difference. Accordingly table 5 collects information on the lowest state in the e and E representations. 
The KB numbers have been left out because they practically duplicate those of Bonn B. The differences in eigenvalues 
between KLS and Bonn B is not large but may have some significance, while the agreement in eigenvectors is complete. 
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Very much the same is true in the comparison of W with the realistic values except for T = 20 and 30, which are not 
among the dominant states of the collective Hamiltonian. The reading of the result is simple: W has discovered and 
dealt with some local problem by altering radically some special matrix elements, but it has scrupulously respected 
the truly important contributions. 

The local problem has several manifestations. In 22 Na the realistic interactions produce a close doublet J = 1, 3, 
instead of the well detached J — 3 observed ground state 0, the 7 band in 24 Mg is too low, and more interesting for 
our purpose, the JT = 20 state in 18 F always too low by about 1 MeV. (It should be noted that the experimental 
counterpart of this level is the second JT = 20 in 18 F, the first is an intruder). 

The solution proposed by the W fit is such that there is no hope to reconcile it with the realistic forces. The 
problem is then to reconcile the latter with the data. We do not know how to do that, and we close this section by 
looking for some clues. 

Though realistic interactions are similar, they also differ in some respects as can be gathered from fig. 3 of ref. Jl9| , 
that shows different calculated spectra of 18 F and 18 0. There are two effects. 



• Overall Dilation. If we refer to eq.(l.E ) the discrepancy can be absorbed by simply changing slightly the oscillator 
frequency, which is in principle fixed by the size of the nucleus. Because of their bad saturation properties this 
choice is not always possible for the realistic potentials. As long as saturation properties have to be treated 
phenomenologically, the overall dilation must be treated as a free parameter. 

• The JT — 01 and 10 matrix elements. The differences due to overall dilation can be eliminated by normalizing 
the interactions to the same a 1 in eq.(2.1). When this is done the most notable discrepancies come from the 
JT = 01 and 10 matrix elements, that are spectroscopically the lightest in the sense that [JT] is smallest. 
As we have seen in the preceeding section they are also the ones most severely affected by renormalization 
uncertainties. 

Whether the local problem in the sd shell can be solved by playing on these different elements is an open question, 
but the analysis in terms of the diagonal representations makes clear that there is no much room for tampering with 
the dominant terms in the interactions. 



APPENDIX A: 



Properties of the f matrices 

There is nothing special about the W matrix that is diagonalized in the E representation except that it is traceless 
by construction, while the / matrix leading to the e representation has a number of non trivial properties. 
Direct properties 

• The / At matrix has twice as many elements as the W JT matrix with J = A: 

if r 7^ s and t^u, W^ stu goes into ujj tsu and toJ ust and the allowed values of T and 7 are the same. 

if r = s or i = u, W^ stu goes into uj^ tsu = uj'^. ust , but T is allowed for (— ) r = —1 only and there 
is no restriction on 7. E.g.: there are 5 possible ways of constructing J = 2 states in the sd shell 
^5/2 > ^3/2 > ^5/2 s i/2> ^5/2^3/2, Si/2^3/2 and there are five JT — 21 states and three JT = 20 states (the first two 
are Pauli forbidden). For the 7 matrices there are 8 possible combinations: the 5 above plus Si/ 2 d 5 / 2 , ^3/2^5/2 
and d 3 / 2 Si/2 which are counted as different. Therefore in W we have a 3 x 3 (r = 20) and a 5 x 5 (r = 21) 
matrix; while in / we have twice an 8 x 8 matrix. Each of these 8x8 matrices consists of two blocks of 3 x 3 
and 5x5 because of the following property. 

• In the e representation we can write 

H = \ E [(fttsu + (-r^; t „J(<Sr 7 A 7 J° + {Prtsu - {-y^frtusWM ] (Al) 

r<t 

S < UT 

in terms of the S and A operators defined in eq.(1.31) of [1]. In the preceding example it means that out of 8 
possible 7 = 2t combinations 3 are of A type and 5 of S type. Note that the matrices are always symmetric 
but the eigenvectors are linear combinations of states that are either of iS or A type. 
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The trace of the / matrix always vanishes, as seen from 



E OrtM = E PrrPu^i-r^ { \ \ \ } = 

•V rt rt ,,m ^ ■> 



7,rt rt 7 [r] 



E^w^« = ( A2 ) 



where we have used (1.1c) and sum rule (A. 19) from [1]. The matrix element W®® tt is always zero because it 
violates Pauli. 

• The variance of the / matrix is a 1 - = -^y-Cvv, where erJV is the variance of the W matrix. The calculation is 
better conducted in the m-scheme 

H(W) =Y,W XX .Z+Z X , 
x=(ij),x = l---D^/2, i<j (A3) 

H(f) — E faa'S a S a ' 

a = (ij) , x = 1 ■ ■ ■ D 2 (ij unrestricted) (A4) 

A typical W XX ' — say, appears as f aa ' = ±-?W xx < in four contributions: aa' = 13 24 and 24 13, 14 23 and 

23 14. Then 



n(2) \ _1 

E w -'; ( A5 ) 



2 

W XT , \ 2 D-l 



Asymmetric factorizations 

Quadrupole renormalizations are mediated by terms of the form k(q' p ■ q' 2huJ + q' 2hw ■ q p ). To understand how 
asymmetric factorization are possible we study the spectra of / x / matrices f(n), whose non zero elements belong to 
the rectangular blocks f xa and f ax , a=l...K,x = K + l...I, I — K = L. 

Specializing eq.( 1.2b| ) to this situation we have 

^ fax u xk = U a kCk ^ fxa.Uak — U x kCk (-^-7) 

x a 

The eigenvector | A; > with eigenvalue can be expanded in terms of unit column vectors \i > (1 in the i-th position, 
zero in the others): 

I ft >= ^u tk \i>= E u a k\a>+ E u xk\x> (A8) 

i=l,I a=l,K x=K+l,I 



By reversing simultaneously the sign of e/. and u x k (&k — ► — &ki u xk — * -"nt, V x) eq.(A7) remains unchanged 
telling us that 

\k > = E u afcl a > -^u xk \x > (A9) 

a x 

is an eigenvector with eigenvalue — e&. Furthermore, from unitarity 

E U ikUik> = E U ak u ak> + E U xkU x k' = <W (A10) 
i a a: 

and taking overlaps 
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< k\k' >= ^UakUak' - Y UxkUxk ' = (All) 
a 

leading to 

"^UakUak' = ^ U x kU x k> = -Skk' (A12) 

a a: 

From these results we may construct the spectrum of f(n). Let's call M = mm(L,K). We have: M positive 
eigenvalues (k — 1, M), M negative ones eg = — ej. (fe = M + k) and 7 — 2M null ones. 

We may gain further insight by presenting the problem as a search for an optimum approximant g a gx to the 
rectangular matrix f ax , defined through 

yX/xa ~ 9a9x) 2 = min . (A13) 
Variation with respect to f x and /„ leads to 



and given that ( A15 ) is invariant under g x — ► ag x g a — > a 1 g a we may request 



and (A14) becomes (A7) by identifying g x = u x k, g a — u a k, e = efe. The factorization produced by the lowest 
(k = 1) eigenstate at e = — |ei| is identical to the one for the highest at e = | e 1 1 and the best available. Exact 
separability is achieved for = for k =/= 1. 

APPENDIX B: 

Baranger and Kumar revisited. 

In two famous papers Baranger and Kumar attempted to derive from a realistic interaction the pairing plus 
quadrupole forces adapted to a space of two major shells pO) and proceeded to do Hartee-Fock-Bogoliubov (HFB) 
calculations in the rare earth region that showed for the first time that it was possible to explain microscopically the 
onset of deformation plj . 

The success of the calculations is probably due in large part to the fact that the model is far more realistic than its 
authors believed. The reason is rather strange and we think it is worth telling. 

Let us start by comparing the traditional pairing plus quadrupole forces as used in |p(i| , pT| 

77p = -G x (P+ + P+ +1 )-(P p + P p+1 ) (Bl) 

x' 

H q = -— {q p + q P +i) ■ {q P + q P +i) (B2) 



with the normalized versions in eqs.( |l.9| ) and 1.10), which we write explicitly by borrowing numbers from table II and 
remembering that hujQ =9 MeV: 

77p = -0.327k.; ( -SL + -p^=j ' ( -5= + -fe= J (B3) 

77, = -0.2167iu, ( * + pL) . f^L + pL) . (B4 ) 
q \JV P M p+ iJ \JV P M p+ iJ 

(B5) 

If we consider first the case of one shell, the operators are the same, and we can relate the coupling constants by 
simply equating. 

For the two shell case, the overlap 
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_L < M. _l qp+1 

V2 N p N p+ i 



qp + qp+1 >=^^±^ = 0.98forp = 3 (B6) 



seems large enough to identify iJ^ and to a good approximation. Of course there is some risk because the norms 
defined in the vector space of the multipolc representation the S r t operators are all treated as unit vectors. In a 
midshell situation, as in 28 Si, the p + 1 - upper - shell is very poorly represented in the wavefunctions. Then the large 
overlap is meanigless. 

With this proviso in mind we equate the traditional and the new forms. Recalling that hui = 40 A -1 / 3 [p2[ , and 
writing the norms in terms of A mp as in eqs.(6,7) we find for one shell 



0-216^ ^ 1 216 ^ = xf, 5/3 0-32^ ^ 19.51 = 2/3 

K* 2 A i/s^V8 2"2 ' fi p A x/3 A V3 " o 



(B7) 



and for two shells, 



0.216(2M ^1 216 _X / _X(, ,- 5 /3. 0-32(2^;) ^ 19.51 2/3 , , 

-^ 2 +AA p 2 +1 " 2 A i/3 A 4/3 2 ~2 A ' n p + Sl p+1 ~ A V*A% 3 U ~ W ^ } 

where we have expressed the averaged norms in terms of A c / p and A c u p , both close enough to the the total number 
of particles at the closure of shell p, A cpi to be identifyed with it in what follows. 

For one shell the only problem comes from the conventional scalings, in A~ 5 ^ 3 and A~ 2 / 3 . It is possible to understand 
their origin: they amount to set A mp = A, which makes sense in comparing strengths in distant regions, but is locally 
wrong. If taken at face value, the A~ 5 ^ 3 local behaviour predicts variations in the moments of inertia of neighbouring 
that are much larger than the observed ones |23|j24| ■ 

In the two shell case the couplings are reduced with respect to the one shell values by a factor of about 

(A mp /A cp f/* = (B9) 

with k = 4 for x' & n d k = 2 forG. If our ide ntifi cations where correct these discrepancies should not exist, and they 
are related to the risk we described after eq. ( p36[). There is no guarantee that we can approximate well the operator 
<?p + Qp+i by q p + q p +i- The compromise x' in ( |B8| ) is too small for the lower shell but also too large for the u ppe r 
one. If the mixing is strong and both shells contribute equally to the wavefunctions, then the large overlap in ( p36] ) 
indicates that the compromise may work but in a nucleus well described by the lower shell alone it makes no sense. 
The thing to do in this case is to restrict the model to one shell. 

By overlapping the q ■ q form with a realistic interaction Baranger and Kumar had obtained x'o = 203 vs our 216 
for one shell, reduced by a factor of about 0.6 for two shells at p = 2, close to the value of 7/8 calculated in eq.(p3S|). 
When faced with this unwanted reduction, instead of blaming the q ■ q form, they declared incorrect their method 
of extraction which - though in a primitive form - is identical to ours and correct to within some details (discussed 
in section II). Then they proceeded to invent another method to obtain the coupling constant, unrelated to any 
interaction (see Note 2 at the end of the section). 

At this point ref. pp| becomes confusing because it is argued that since the quadrupole force cannot be extracted 
from the realistic interaction its origin must be something else that may not be quadrupole at all. What we are showing 
is that the "something else" is simply the normalized quadrupole force. 

The extraordinary thing is that Baranger and Kumar had found it! The reason they did not see what they had 
found is that they reasoned in terms of the A~ 5 ^ 3 scaling, in spite of having given the correct argument about what 
the scaling must be! 



Now observe carefully eqs.(B7) and ( JB8| ) to discover the miracle that repairs the dammage. The only sensible way 
to define x'o f° r one shell is to "equate" A mp « A, while for two shells we must take A cp « A and now x'o is identical 
to 216 in both cases! 

There is no contradiction between 

X 1 (1 shell)=/= x' '(% shells) at a given nucleus and 

Xo(l shell)= x'o(<? shells), both calculated in different nuclei. 

This is the subttle argument missing in [ pp[ It is only thanks to universal scaling that it can be made. 

Everything seems to happen as if one flaw of the model - incorrect scaling - corrected the other: the space dependence 
of the coupling constants. A better interpretation though, is that the model should be restricted either to one shell 
in the vicinity of A mp , or to two shells in the vicinity of A cp . 
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Now we note that the rare earth region studied in |2l|] is approximately centered at the oscillator closures Z=70 
and N=112, which meets this restriction, and that the parameters used were (in Mev): 

Xo ~ 280, Gqtt = 27 Gq v = 22, 



somewhat larger than those in eq.(B8) but quite consistent with the renormalized values of section III. ( The need 
to use different Go values for neutrons and protons is readily explained by ( p38| ). It is a mild manifestation of space 
dependence) . 

We can draw two conclusions: 

i) the conventional pairing plus quadrupole model in two shells is far more consistent with the realistic interactions 
than its authors believed. 

it) it can be made truly realistic by using the normalized operators, by including a pair of other terms (octupole 
and hexadecapole), and by examining more closely the monopole contribution. Very little of the basic simplicity of 
the original will be lost in this improved version. 

Note 1. On space dependence. 

Examine what happens in larger spaces when eq.(fB§|) is generalized to M major shells. It is elementary to prove 
that x' = 0(A~ 1/3 M- 4 ), and the overlap tends to and x' = 0(A-^ 3 M- 4 ). It means that to simulate - never 

mind how remotely - the behaviour of its realistic counterpart, the conventional q ■ q force must be affected by a 
vanishingly small x' ■ 

Note 2. On the second method of extraction in fi2l\]. To replace direct extraction of x' from the realistic interaction 
Baranger and Kumar proposed a method based on the idea that the average energy of a nucleus is independent of its 
shape. It makes no reference to the force nor to any empirical datum. It leads to: 



X 



(eqs. (72,73) and paragraph followin g th em in pOj j, but here we have introduced the norms instead of their asymptotic 
values). This is the same x' °f cq.(|B8|), except that a factor 4 x 0.216 has been set to unity. The identity in form 
suggests that the method is a valid dimensional analysis, but the similarity in the numbers is likely to be accidental. 



TABLE IV. Comparison of some matrix elements in the sd shell (5 
indicates large W-realistic discrepancies. 



d 



5/2, 



d 



3/2; 



1 = S1/2) for different interactions. 



rstu 


JT 


KLS 


KB 


Bonn B 


W 


5553 


10 


3.03 


3.17 


3.31 


2.54 




21 


-0.52 


-0.40 


-0.22 


-0.28 




30 


1.21 


1.87 


1.89 


2.22 




41 


-1.24 


-1.36 


-1.28 


-1.24 


5533 


01 


-4.17 


-3.79 


-3.41 


-3.19 




10 


1.45 


1.62 


1.29 


0.72 




21 


-0.89 


-0.90 


-0.89 


-1.62 




30 


0.13 


0.50 


0.56 


1.89 <- 


5153 


20 


-1.10 


-1.44 


-1.33 


-0.1 <- 



TABLE V. Lowest energies Er for L = 01, 10, 20 and 30 and e 7 for all 7, and overlaps of the wavefunctions for KLS, BonnB 
and W. 



01 


10 


20 


30 




10 


11 


20 


21 


30 


31 


40 


41 


-5.42 


-5.43 


-2.68 


-2.15 


KLS 


-2.18 


2.38 


-2.90 


-0.71 


-0.82 


0.44 


-1.61 


0.40 


-5.48 


-6.24 


-2.91 


-2.66 


BonnB 


-1.55 


2.64 


-3.32 


-0.97 


-0.83 


0.46 


-1.39 


0.52 


-5.69 


-5.90 


-0.95 


-2.44 


W 


-2.16 


3.08 


-3.18 


-0.70 


-0.94 


0.54 


-1.60 


0.51 


1.00 


1.00 


1.00 


1.00 


<KLS|BonnB> 


.99 


1.00 


1.00 


.99 


.98 


.99 


1.00 


1.00 


1.00 


0.98 


0.55 


0.82 


<KLS| W> 


.95 


.99 


1.00 


.98 


.98 


.92 


1.00 


1.00 
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FIG. 1. ^-eigenvalue density for the KLS interaction in the pf+sdg major shells hco = 9. Each eigenvalue has multiplicity 
[r]. The largest ones are shown by arrows. 
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FIG. 2. e-eigenvalue density for the KLS interaction in the pf+sdg major shells. Each eigenvalue has multiplicity [7]. The 
largest ones are shown by arrows. 
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FIG. 3. iS-eigenvalue density for the KLS interaction in the pf+sdg major shells hco — 9, after removal of the five largest 
multipole contributions. Each eigenvalue has multiplicity [F]. The largest ones are shown by arrows. 
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FIG. 4. E-eigenvalue density for the KLS interaction in the pf+sdg major shells hu> = 9, after removal of multipole contri- 
butions with \e\ > 1.3. Each eigenvalue has multiplicity [r]. The largest ones are shown by arrows. 
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